Existence of eigenvalues of a linear operator pencil in 
a curved waveguide — localized shelf waves on a 

curved coast 

E R Johnson 
Department of Mathematics, University College London 
Gower Street, London WCIE 6BT, U. K. 
email erj@math.ucLac.uk 

Michael Levitin 

Department of Mathematics, Heriot-Watt University 
Riccarton, Edinburgh EH14 4AS, U. K. 
email M.Levitin@ma.hw.ac.uk 

Leonid Parnovski 

Department of Mathematics, University College London 
Gower Street, London WCIE 6BT, U. K. 
email Leonid@math.ucl.ac.uk 

September 2004* 



Abstract 

The study of the possibility of existence of the non-propagating, trapped 
continental shelf waves (CSWs)along curved coasts reduces mathemati- 
cally to a spectral problem for a self-adjoint operator pencil in a curved 
strip. Using the methods developed in the setting of the waveguide 
trapped mode problem, we show that such CSWs exist for a wide class of 
coast curvature and depth profiles. 
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1 Introduction 

Measurements of velocity fields along the coasts of oceans throughout the world 
show that much of the fluid energy is contained in motions with periods of a few 
days or longer. The comparison of measurements at different places along the 
same coast show that in general these low-frequency disturbances propagate 
along coasts with shallow water to the right in the northern hemisphere and 
to the left in the southern hemisphere. These waves have come to be known 
as continental shelf waves (CSWs). The purpose of the present paper is to 
demonstrate, using the most straightforward model possible, the possibility of 
the non-propagating, trapped CSWs along curved coasts. The existence of such 
non-propagating modes would be significant as they would tend to be forced by 
atmospheric weather systems, which have similar periods of a few days, similar 
horizontal extent, and a reasonably broad spectrum in space and time. Areas 
where such modes were trapped would thus appear to be likely to show higher 
than normal energy in the low frequency horizontal velocity field. 

The simplest models for CSWs take the coastal oceans to be inviscid and of 
constant density. Both these assumptions might be expected to fail in various 
regions such as when strong currents pass sharp capes or when the coastal 
flow is strongly stratified. However for small amplitude CSWs in quiescent flow 
along smooth coasts viscous separation is negligible. Similarly most disturbance 
energy is concentrated in the modes with the least vertical structure, which are 
well described by the constant density model [LBM]. The governing equations 
are then simply the rotating incompressible Euler equations. Further, coastal 
flows are shallow in the sense that the ratio of depth to typical horizontal scale 
is small. Expanding the rotating incompressible Euler equations in powers of 
this ratio and retaining only the leading order terms gives the rotating shallow 
water equations [PeJ. 

du ~ 

(1.1) — -I- u • grad u — 2i7k x u = —g grad H, 

r)J-f ~ 

(1.2) _ + div [(F + /7)u] = 0. 

Here x and y are horizontal coordinates in a frame fixed to the rotating Earth, 
k is a vertical unit vector, u{x,y,t) is the horizontal velocity, Q, is the (locally 
constant) vertical component of the Earth's rotation, g is gravitational acceler- 
ation, H{x,y,t) is the vertical displacement of the free surface and H{x,y) is 
the local undisturbed fluid depth. 

System ()l.lj) . ()1.2j) admits waves of two types, denoted Class 1 and Class 

2 by [LaJ. Class 1 are fast high-frequency waves, the rotation modified form 
of the usual free surface water waves, although here present only as long, non- 
dispersive waves with speeds of order \/gH. Class 2 waves are slower, low- 
frequency waves that vanish in the absence of depth change or in the absence 
of rotation. It is the Class 2 waves that give CSWs. They have little signature in 
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the vertical height field H{x, y, t) and are observed through their associated hor- 
izontal velocity fields | Ha| . The Class 1 waves can be removed from (|l.l|) . p.2|) 
by considering the 'rigid-lid' limit, where the external Rossby radius ^gH 
(which gives the relaxation distance of the free surface) is large compared to 
the horizontal scale of the motion. This is perhaps the most accurate of the 
approximations noted here, causing the time-dependent term to vanish from 
()1.2jl and the right sides of (fl.ljl to become simple pressure gradients. 

For small amplitude waves the nonlinear terms in ()l.lj) . ()1.2j) are negligible 
and cross-differentiating gives 

(1.3) -^ + 2f7divu = 0, 

(1.4) div(iJu) = 0. 

where C = ~ f| is the vertical component of relative vorticity. Equation 
4]) is satisfied by introducing the volume flux streamfunction defined though 



allowing p.3jl to be written as the single equation 

(1.6) div ( ^gi'ad ) + • grad ^ x grad {XjH) = 0. 

\li Ot J 

Equation ()l.bj) is generally described as the topographic Rossby-wave equation 
or the equation for barotropic CSWs. Many solutions have been presented for 
straight coasts, where the coast lies along y = (say) and the depth H \s a 
function of y alone (described as rectilinear topography here) [LBMJ. These 
have shown excellent agreement with observations of CSWs, as in [HaJ. There 
has been far less discussion of non-rectilinear geometries, where either the coast 
or the depth profile or both are not functions of a single coordinate. Yet inter- 
esting results appear. [StHull l5tHu2| present extensive numerical integrations 
of a low-order spectral model of a rectangular lake with idealized topography. 
For their chosen depth profiles normal modes can be divided into two types; 
basin-wide modes which extend throughout the lake and localized bay modes. 
These bay modes correspond to the high-frequency modes found in a finite- 
element model of the Lake of Lugarno by [TrJ and observed by [StHu SaTrZa| . 
\ Jo2. StJoll l5tJo2 j give a variational formulation and describe simplified quasi- 
analytical models that admit localized trapped bay modes. However the ge- 
ometry changes in these models are large, with the sloping lower boundary 
terminating abruptly where it strikes a coastal wall. Further [JolJ notes that 
(|1.6|1 is invariant under conformal mappings and so any geometry that can be 
mapped conformally to a rectilinear shelf cannot support trapped modes. The 
question thus arises as to whether it is only for the most extreme topographic 
changes that shelf waves can be trapped or whether trapping can occur on 
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smoothly varying shelves. The purpose of the present paper is to show that this 
is not so: trapped modes can exist on smoothly curving coasts. 

The geometry considered here is that of a shelf of finite width lying along 
an impermeable coast. Thus sufficiently far from the coast the undisturbed 
fluid depth becomes the constant depth of the open ocean. It is shown in [Jo3j 
that at the shelf-ocean boundary of finite-width rectilinear shelves the tangential 
velocity component u vanishes for waves sufficiently long compared to the shelf 
width. The wavelength of long propagating disturbances is proportional to their 
frequency which is in turn proportional to the slope of the shelf. Thus it appears 
that for sufficiently weakly sloping shelves the tangential velocity component, 
i.e. the normal derivative of the streamfunction, at the shelf-ocean boundary 
can be made arbitrarily small. Here this will be taken as also giving a close 
approximation to the boundary condition at the shelf-ocean boundary when this 
boundary is no longer straight. The unapproximated boundary condition is that 
the streamfunction and its normal derivative are continuous across the boundary 
where they match to the decaying solution of Laplace's equation (to which ()1.6jl 
reduces in regions of constant depth). This gives a linear integral condition along 
the boundary. The unapproximated problem will not be pursued further here. 
The boundary condition at the coast is simply one of impermeability and thus on 
both rectilinear and curving coasts is simply that the streamfunction vanishes. 
Now consider flows of the form 

(1.7) ip{x, y, t) = Re{<I>(rE, y) exp(-2iu;i7t}, 

so <^(x,y) gives the spatial structure of the flow and cu its non-dimensional 
frequency. Then $ satisfies 

(1.8) + gi-ad • g^ad ^ + • (^grad ^> x grad (^^^ ^ = 0, 

(1.9) $ = at the coast, 

(1.10) n • grad <I> = at the shelf-ocean boundary, 

where n is normal to the shelf-ocean boundary. 

Mathematically, we are going to study the existence of trapped modes (i.e., 
the eigenvalues either embedded into the essential spectrum or lying in the 
gap of the essential spectrum) for the problem p.8j) - p.l0j) in a curved strip. 
Similar problems for the Laplace operator have been extensively studied in the 
literature — either in a curved strip, or in a straight strip with an obstacle, 
or in a strip with compactly perturbed boundary. In the case of Laplacians 
with Dirichlet boundary conditions these problems are usually called 'quantum 
waveguides'; the Neumann case is usually referred to as 'acoustic waveguides'. 
The important result concerning quantum waveguides was established in [ExSe^, 
|DuEx| : in the curved waveguides there always exist a trapped mode. Later this 
result was extended to more general settings; in particular, in |KrKrjj it was 
shown that in the case of mixed boundary conditions (i.e., Dirichlet conditions 
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on one side of the strip and Neumann conditions on the other side) trapped 
modes exist if the strip is curved 'in the direction of the Dirichlet boundary'. 

The case of quantum waveguides is more complicated because any eventual 
eigenvalues are embedded into the essential spectrum and are, therefore, highly 
unstable. Therefore, it is believed that in general the existence of trapped modes 
in this case is due to some sort of the symmetry of the problem (see |EvLeVa| . 



In the present paper we use the approach similar to the one used in [DuEx| 
and [EvLeVaJ : however, we have to modify this approach substantially due to 
the fact that we are working with a spectral problem for an operator pencil 
rather than that for an ordinary operator. 

The rest of the paper is organized in the following way. In Section 2, we 
discuss the rigorous mathematical statement of the problem; in Section 3, we 
study the essential spectrum, and in Section 4, state and prove the main result 
on the existence of a discrete spectrum (Theorem 4.1). In particular we show 
that a trapped mode always exists of all of the following conditions are satisfied: 
(a) the depth profile H does not depend upon the longitudinal coordinate; (b) 
the channel is curved in the direction of the Dirichlet boundary; (c) the curvature 
is sufficiently small. 

Similar results can be obtained in a straight strip if the depth profile H 
depends nontrivially upon the longitudinal coordinate; we however do not discuss 
this problem here. 

2 Mathematical statement of the problem 

2.1 Geometry 

The original geometry is a straight planar strip of width 5: 



Deformed geometry G is assumed to be a curved planar strip of con- 
stant width 6. To describe it precisely, we introduce the curve F = {{x = 
■^{0^ y = ^(0)}) G I^. where ^ is a natural arc-length parameter, i.e., 
X'iO^ + Y'iO^ = 1 . By 7(0 = X"{Oy'iO - X'iOy'iO we denote a 
(signed) curvature of T (see FigureH]). Note that |7(0P = X"{0^ + Y"{0^. 

We additionally assume 

(2.1) supp7 ^ [— -R, i?] for some i?>0, 



[D^. [AsPaVajV 



Go = {{x,y):xeR, ye (0,6)}. 



and set 



(2.2) 



K 



,+ 



sup 7(0 , 



K 
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Y(s)<0 



Figure 1: Domain G and curvilinear coordinates ^, 77. The solid line denotes 
the boundary diG with the Dirichlet boundary condition, and the dotted line 
— the boundary (?2G with the Neumann boundary condition. 



We shall assume throughout the paper the smoothness condition 

(2.3) 7GC~(M), 

which can be obviously softened. 

Now we can introduce, in a neighbourhood of V, the curvilinear coordinates 

(C,??) as 

(2.4) X = x{0 - vY'iO , y = y{0 + vx'iO , 

(where ?? is a distance from a point {x,y) to F) and describe the deformed strip 
G in these coordinates as 

(2.5) G = G^ = {iC,v) ■■ CeR, r/G (0,<5)}. 

Remark 2.1. As sets of points, G^ = Gq for any 7, but the metrics are different, 
see below. We shall often omit index 7 if the metric is obvious from the context. 

To avoid self-intersections, we must restrict the width of the strip by natural 
conditions 

(2.6) < AS-'^ , A = const G [0, 1) . 

Finally, it is an easy computation to show that the Euclidean metric in the 
curvilinear coordinates has a form dx'^ + dy^ = gd^"^ + drf', where 

5(e,r?) = (l + r?7(0)'- 
Later on, we shall widely use the notation 

(2.7) p(e,r/) = (5(C,77))V2 = i + ^^(^). 

Note that in all the volume integrals, 

dG^ = p(C, V) di dr? = (1 + 777(0) d?? = p{i, v,) dGo • 
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2.2 Governing equations 

For a given positive continuously differentiable function H{^,r]) (describing a 
depth profile), we are looking for a function $(^, rj) satisfying (|1.8j) with spectral 
parameter uj. 

By substituting 

(2.8) /?(^,r/) :=ln//(^,77), 

and using explicit expressions for differential operators in curvilinear coordinates, 
we can re-write (|1.8|) as 

(2.9) 

1 d"^^ ^ / 1 dp ^ 1 dl3\ d<^ ^ /dl3 1 dp\ d<i>\ 



p2 Q^2 Qj^2 p2 J \dr] p Or] J df] J 

p\d^ drj drj dS, 

Remark 2.2. When deducing equation ()2.9jl . we have cancelled, on both sides, 

a common positive factor h(^,r]) := — = e~^^^'''\ However, we have to 

use this factor when considering corresponding variational equations, in order to 
keep the resulting forms symmetric. This leads to a special choice of weighted 
Hilbert spaces below. 

Further on, we only consider the case of a longitudinally uniform monotone 
depth profile 

(2.10) p(^^r,)^P{rj), P'{v)>0, 

in which case the equation ()2.9jl simplifies to 



(2.11) 

with P' = f. 

drj 



^ , 1 92^) a2^> ^ 1 dp d<^ ^ /, 1 dp\ 

p2 Q^2 Q^2 p3 y p dr] J dr] 

p 5^ ' 



2.3 Boundary conditions 

Let diG = {(C,0) : C e and = {(^,<5) : ^ e M} denote the lower 
and the upper boundary of the strip G, respectively. Boundary conditions ()1.9j) . 
(|1.10jl then become 

drj 

Remark 2.3. If the flow is confined to a channel then the Dirichlet boundary 
condition (|1.9|1 applies on both channel walls. This leads to a mathematically 
different problem which we do not consider in this paper. 



(2-12) '^\a,G 



= 0. 

d2G 



Existence of eigenvalues of an operator pencil in a curved waveguide 



Johnson/Levitin/Parnovski — September 2004 



Page 8 



2.4 Function spaces and rigorous operator statement 

We want to discuss the function spaces in which everything acts. Let us denote 
by L2{G;h) the Hilbert space of functions (j) : G —>■ C which are square- 

integrable on G with the weight h{r]) = — = exp(— /3(r/)): 

H 



2 

L2{G;h) 



m,V)\^hi7j)dG 



G 



The corresponding inner product will be denoted ■) L2{G-h)- Similarly we can 
define the space L2{F] h) for an arbitrary open subset F of G. 
Let us formally introduce the operators 



1 a2^> 1 dp 



+ 



p2 Q^2 Q^2 p3 



1 dp\ 



p dr] J dr] 



and 



p 5^ 



(The dependence on 7 is of course via p, see ()2.7jl .) Then ()2.11jl can be formally 
re-written as 

(2.13) uC^^ = M^^, 
or via an operator pencil 

(2.14) Aj = Aj{uj) = iuCj - 
as 

(2.15) A^{io)<^> = 

The domain of the pencil Aj in the L2-sense is naturally defined as 

(2.16) Dom(^^) = e H^{G) , $ satisfies (|ZT^ T , 

where denotes a standard Sobolev space. 

On the domain ()2.16jl . Ai^ is symmetric, and is symmetric and positive 
in the sense of the scalar product (•, •) L^(^G-h)' ^'^'^ 



{G;h) 



+ p 



drj 



Later on, we shall use a weak (or variational) form of (|2.15|1 . and shall 
require some other function spaces described below. Let F C G, and suppose 
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its boundary is decomposed into two disjoint parts: dF = d\F U 32-F. We 
introduce the space 

c^{F■^lF) = {4>ciC^{F)■.^^I^r^^lF = 

and there exists r > such that (/){£,, r]) = for {£,,7]) G F, \^\ > r} , 

consisting of smooth functions with compact support vanishing near diF. By 
Hq{F; diF; h) we denote its closure with respect to the scalar product 

i^^'^)Hl,{F,d,F;h) = i^^^) L2{F;h) + <t^.&^^ i') L^iF.h) ■ 

In what follows we shall study the operators and the pencil 

from a variational point of view. The details are given in the next Section; here 
we note only that from now we understand the expression (^C^^', ^)^2(G^,/i) 
as the quadratic form for the operator with the quadratic form domain 

Hl{G-diG-h). 

The main purpose of this paper is to study the spectral properties of the 
operator pencil A^. We recall the following definitions. 

A number G C is said to belong to the spectrum of A^ (denoted 
spec(^^)) \^ A~f{io) is not invertible. 

It is easily seen that in our case the spectrum of A-y is real. 

We say that w € M belongs to the essential spectrum of the operator pencil 
A^ (denoted to G speCgs5(^-y)) if for this lu the operator ^^(w) is non-Fredholm. 

We say that lu G C belongs to the point spectrum of the operator pencil 
Aj (denoted lu G speCpt(w4^)), or, in other words, say that lu is an eigenvalue, 
if for this lu there exists a non-trivial solution G Dom(^^) of the problem 
A^iuj)"^ = 0. 

It is known that the essential spectrum is a closed subset of M, and that any 
point of the spectrum outside the essential spectrum is an isolated eigenvalue 
of finite multiplicity. The set of all such points is called the discrete spectrum, 
and will be denoted speCcjis(^-y) There may, however, exist the points of the 
spectrum which belong both to the essential spectrum and the point spectrum. 

Our main result (Theorem 14.11 below) establishes some conditions on the 
curvature 7 of the waveguide which guarantee the existence of eigenvalues of 

It is more convenient to deal with problems of this type variationally, and 
we start the next Section with an abstract variational scheme suitable for self- 
adjoint pencils with non-empty essential spectrum. 
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3 Essential spectrum 

3.1 Variational principle for the essential spectrum 
Definition 3.1. We set, for j G N, 

(3.1) fJ'y-i = sup mi , ^ ^ — — ' — 

dim U =j 

As {C^"$,^)L2{G-y,h) is positive, the right-hand sides of ()3.1|1 are well de- 
fined, though the numbers fi^j may a priori be finite or infinite. 

Obviously, for any fixed curvature profile 7 the numbers fiyj form a non- 
increasing sequence: 

/^7,l > ^7,2 > • • • > > /^7,j + l > • • • • 

Definition 3.2. Denote 



(3.2) fj.-y = lim /i^ 



3- 



For general self-adjoint operator pencils the analog of ()3.2j) may be finite or 
equal to ±00; as we shall see below, in our case is finite 

The following result is a modification, to the case of an abstract self-adjoint 
linear pencil, of the general variational principle for a self-adjoint operator with 
an essential spectrum, see || Dal Prop. 4.5.2]. 

Proposition 3.3. Either 

(i) fi^ > —00, and then supspeCe55(.A-^) = n-y, 

or 

(a) fi^ = —00, and then sY>eCgs^{A^i) = 0. 

Moreover, if fJ'-yj > fJ-y, then jj^^j G spec^i^{Ay). 

Proposition 13.31 ensures that we can use the variational principle p.lj) in 
order to find the eigenvalues of the pencil lying above the supremum {j.^ of 
the essential spectrum. 



3.2 Essential spectrum for the straight strip 

The spectral analysis in the case of a straight strip (7 = 0) is rather straight- 
forward as the problem admits in this case the separation of variables. 

Let us seek the solutions of ()2.11jl . ()2.12j) in the case of a straight strip 
(7 = 0, and so p = 1) in the form 

(3.3) «>(e,7?) = </'(r/)exp(iae); 
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it is sufficient to consider only real values of a. 

After separation of variables, (|2.11|) . (|2.12|1 are written, for each a, as a 
one-dimensional transversal spectral problem 

(3.4) u;(-(/>" + /3V + aV) = a/3>, (/>(0) = = , 

Alternatively, introduce operators 

(q, : 1-^ -(/)" + + a^^ , ma', aP'(j) , 

and a pencil 

(again understood in an L2((0,5);/i) sense with the domain defined similarly 
to ()2.1&j) ). and consider a one-dimensional operator pencil spectral problem 

aa{oj)(i) = 0. 

For a fixed value of a, the one-dimensional linear operator pencil (|3.4|) 
has the essential spectrum {0} and a discrete spectrum spec(aa); note that 
spec(a_o) = — spec(oo). Denote, for a > 0, the top of the spectrum of this 
transversal problem by Ua = supspec(oQ,). 

Lemma 3.4. Let a > 0. Then, under condition ()2.10j) . 

(i) spec(oQ) C [0, +co); 

(ii) < LOa < +oo; 

(Hi) ujce +0 as a ^ oo. 
Proof. By the variational principle analogous to Proposition [^fi), 

Wo = ^ sup Ja{(p) , 

where we set 
(3.5) 

5 



J (-0) _ ("^a<^''^)i2{(0,5),h) _ 



-(l)"{ri) + P'{r])(l)'{ri) + a'^(f){ri))(p{r])h{r]) di] 





After integrating by parts using h{rj) = e~^^^\ and inverting the quotient, 
we get 

(3.6) j„(0)= + 

\ a 
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where we denote 



and 



I /3'(r?)e-'^('')|0(ry)|2d7? 



-^('')|(/)'(r/)|2d?7 



/3'(r/)e-'^('')|</.(r?)pd?7 



The statements (ii) and (iii) of the Lemma now follow immediately from the 
estimates 

inf e-^^''^ 
j(i)(<^)>^^ 



and 



sup (/3'(r?)e-'^('^)^ 

'76(0,5) 

inf e-^^'') 



4(^2 sup (/3'(r/)e-^('')) ' 

where the latter inequality uses the variational principle and the fact that the 
principle eigenvalue of the mixed Dirichlet-Neumann problem spectral problem 

for the operator —-r^ on the interval (0,6) is equal to —pr. The statement (i) 

dr/^ 4o^ 

follows from the positivity of the right-hand side of p.6jl . 



□ 



We are now able to find the essential spectrum of the problem in a straight 
strip. 

Lemma 3.5. Assume that conditions ()2.1Qjl liold. Then 
(3.7) speCess(A) = , 

where 



/5'(^)e-/3W|0(^)|2dry 



(3.8) = sup 







> 0, 



S 5 
e-^('')|0'(7?)|2dr/-y" e-^('')|0(r?)|2dr? 
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Proof. It is standard that 

speCe5s(A) = IJ spec(aQ,) . 

am. 

Thus, by Lemma IT4l and with account of the anti-symmetry of the spectrum 
of Oq, with respect to a and its positivity for a > 0, we have 

SUpspeCgs5(^o) = SUpWa = sup sup Jai^i) ■ 

a>0 a>0</,g|/i((o,5),0,/i) 

By maximizing first with respect to a, we obtain, from p.6jl . 

Ja{(l>) < 

with the maximizer 



«*(</>) = 




Maximizing now with respect to gives supspeCgs5(^o) = with 
given by ()3.8jl . Finally, the statement follows from the fact that cua depends 
continuously on a > and thus takes all the values in (0, ri*] by Lemma lT4l □ 



3.3 Essential spectrum for a curved strip 

It is now a standard procedure to show that under our conditions the essential 
spectrum of the problem in a curved strip coincides with the essential spectrum 
of the problem in a straight strip given by Lemma ITSl Namely, we have the 
following 

Lemma 3.6. Let us assume the conditions (fZT|l . (fZ^ and (fZTn|l iiold. Then 

speCess(-^7) = speCess(.Ao) = 

with given by ()3.8j) . 

The proof is based on the fact that any solution of the problem ()2.11jl . 
(fZT^ with 7^0 (and thus p ^ 1) should coincide in G n {|^| > R > 
max(| inf supp7|, | supsupp7|)} with a solution of the same problem for 7 = 0. 
An analogous result has been proved in a number of similar situations elsewhere, 
see e.g. |ExSel lEvLeVal fPaPaj . so we omit the details of the proof. We briefly 
note that the inclusion speCes5(^^) C speCe5s(.Ao) is proved using the separa- 
tion of variables as above and a construction of appropriate Weyl's sequences, 
and in order to prove the inclusion speCes5(^^) 5 speCg5s(^o) one can use 
the Dirichlet-Neumann bracketing and the discreteness of the spectrum of the 
problem (EHU), (|TT^ considered in Gn{|^| < R} with additional Dirichlet or 
Neumann boundary conditions imposed on the "cuts" = ±R}. 
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4 Main result 

Our main result consists in stating some sufficient conditions on the depth 
profile (3{ri) and the curvature profile 7(^) which guarantee the existence of an 
eigenvalue of the pencil Aj lying outside the essential spectrum. 

Theorem 4.1. Assume, as before, that the condition ()2.10j) holds. Assume 
additionally that 

(4.1) I3"{v)<0 fori] e (0,5). 

Then there exists a constant Cp > 0, which depends only on the depth profile 
P, such that spec^is{A.y) / whenever j satisfies conditions ()2.1j) . ()2.3j) and 

(4.2) 1 7(e)de yc^l 7(e)'de 

We give an explicit expression for Cp below, see ()4.15jl . 
An integral sufficient condition (|4.2|1 may be replaced by a pointwise, al- 
though more restrictive, condition: 

Corollary 4.2. Assume that the conditions (jm?]) and (|4rT|l hold. Then there 

C 

exists a constant Cj^^R = which depends only on the depth profile (3 and 

2R 

a given R > such that specj;5(.A^) ^ whenever 7^0 satisfies conditions 
(HU), JZD and 

(4.3) 0<7(e)<c^,fl for\C\<R. 

We prove Theorem 14.11 using a number of simple lemmas, the central of 
which is the following 

Lemma 4.3. Suppose there exists a function \I' G Hl{G^,h) such that 

(4.4) (WW.,.)^^^^ 

Then there exists w > il* which belongs to specj,5(.A^). 

Lemma 14.31 is just a re-statement of the variational principle of Proposi- 
tion 13.31 The main difficulty in its application is of course the choice of an 
appropriate test function ^ . However such choice becomes much easier if we 
use the following modification of this Lemma which allows us to consider test 
functions which are not necessarily square integrable on G^. 

Denote, for brevity, G!^ = n {|^| < r}. 

Lemma 4.4. Suppose there exists a function ^ and a constant D such that, 
for any r > R, we have ^ G H^{Gl^, h) and 

(4.5) (Al^vl/, ^)L,{Gr,h) - ^*{C^^, ^)L2(G:;,h) >D>Q. 
Then there exists w > il* which belongs to specj,5(^^). 



Existence of eigenvalues of an operator pencil in a curved waveguide 



Johnson/Levitin/Parnovski — September 2004 



Page 15 



The proof of Lemma 14.41 uses the construction of an appropriate cut-off 
function x(0 such that ^' = satisfies the conditions of Lemma 14.31 cf. 
|DaPal Prop. 1]. 

We now proceed as follows. 

Let be a minimizer in p.8jl . and set 

where 



(4.6) Q, = a*(0*) 



It is important to note that is in fact an "eigenfunction" of the essential 
spectrum of corresponding to its highest positive point il* and that (j)^ is an 
eigenfunction of p.4j) with a = a, (i.e. of the pencil Oq. ) again corresponding 
to the eigenvalue Q^, and so 

(4.7) = + (a^ - A,a./3')<^* , <^.(0) = (5) = 

with A* := (cf. (HH)). 

For future use, we summarize the relations obtained so far: 

Co^ = (-0';(r?) +/3'(r?)0;(r/) +a.</.,(7?))e*-« = ([„.</..) e*"-« , 
Mo^ = a.(3'{r,)Mv)e'"'^ = (m,.0.) e^°-« , 

p{i,n) = 1 + 7(6^. 

(with ' denoting differentiation with respect to ry). 
It is important to note that for any r > 

(A^0^,^)L2(Gg,fe) _ {ma,4'*,(t^*)L2{(0,S),h) 
(£0^, ^)L2{Gl,h) 0*)L2({O,<5),h) 

and, as explicit formulae above show, 
(4.8) 

r 5 



17* > 



-r 

U j /3'(r?)e-^(''V*(^)l'dry>0. 



2rQ, 
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We want to show that under conditions of Theorem 14. H and with the choice 
of ^' as above, the inequality (|4.5j) holds for any r > R. 
In view of ()4.8jl . it is enough to show that 

is negative for r > R. 

Explicit substitution gives, after taking into account the formula 

1 ^^'(^'^).d^ = o 



(due to ()2.1jl . with account of ()2.7jl ). the following expression 

r S r 5 

-r -r 

This, in turn, can be re-written, using the obvious identity 



as 



(4.9) 



r S 

D, = J 7(0 / Ve-^^''^ mv)\' - al\Mv)\') dr?d^ 
-r 

-r 



We shall deal with the two terms in ()4.9j) separately. 
The first one is more difficult. As ()4.6jl yields explicitly 



5 



2 

a. — 



<5 
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we get 
(4.10) 

<5 

h:= y"?7e-^('') {WM?-al\<PM?) dry 



= ^ X I / 7?e-'^('^)|0:(r/)|2dr?. / e''^^'^) |<A.(r/)|2 dr? 



0/ 

We want to show that the term in brackets is negative under some reasonable 
assumptions: 

Lemma 4.5. Assume that the conditions of theorem \4.1\ hold. Then Ii < 0. 

The proof of Lemma |431 uses the following simple fact:^ 

Lemma 4.6. Let (a, b) C (0, +oo) be a finite interval, and let a function 
g : {a,b) —>■ R be non-negative and non-increasing. Then 

xg{x)f{x) dx I • I / f{x) dx I - I / g{x)f{x) dx • / xf{x) dx I < 



for any non-negative function / : (a, 6) — > M. 
Proof of Lemma We have 

•■"'■■»)•(/'-)-(/"'-) •(/■'- 

b h b b 

xg{x)f{x)f{y) dxdy - j j 9{x)f{x)yf(y) dxdy 

a a a a 

by b b 

{x - y)f{x)f{y)g{x) dxdy + j y)f{x)f{y)g{x) dxdy 

a y 



^we are grateful to Daniel Elton for a useful suggestion helping to prove this Lemma. 
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Interchanging the variables x and y in the last integral, we obtain that the whole 
expression is equal to 

h y 

f{x)f{y){g{x) - g{y))^ dxdy , 

o- (I non-positive non-negative 

and is therefore non-positive. □ 

We can now proceed with evaluating Ii. 

Proof of Lemma [T51 We act by doing a lot of integrations by part. We shall 
also use (|47|l . 

We have (all integrals are over [0,6] and with respect to i]) 



J ■ {e-W* + {Oil - A*a.P')e-^V<P*) ■ 



Further, 



*e 



thus producing 

-/3U'|2 



(4.11) negative constant 



Also, 



(4.12) =- j e-^* (-/?>'* + + - A,a,(3'c^,) 

= - j e-^cl)l {al - A,a./3') . 
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Substituting (|4nHl and ^AT^ into ^AJU^ . and simplifying, we get 
(4.13) 

/i-e-^('')|(/>.(r/)|2= f 



positive integral 



-/3U'|2 



?7e 



-/3U |2 



\ negative as /3'>0 negative constant / 



+ve constant ^ 



'.-/5u |2. / e-/3u |2 



/3'e-^|(/.,|2. /7?e-^|0,|2 



non-positive by Lemma f4.6l vtfith g=f3' , f=e as g'=/3"<0 



Thus /i < 0. 



□ 



Let us now return to ()4.9jl and deal with the second term in the right-hand 
side. We have, with account of and (fZ^ . 



1 + r/7(0 
and so 

where 



< < 1 



l-A 



if 7(e) >0 I 



r 5 



-r 



(4.14) 



Jo := max] 



{l,-^}al [ ryV^('')|<A*(r/)|2dry>0. 



1 -^^ 

Thus, as 7(e) vanishes for \^\ > R, we have 



D^ = IiJ 7(e)de + /2 y 7(0'dC = (-/i) y 7(0'dC - J 7(0dc) , 



where 



(4.15) 







r^e-^^") (|</.:(7?)|2-a2|<^,(,?)|2) d,? 
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is a positive constant. 

As soon as ()4.2j) holds, is negative, and so (|4.5j) holds. This proves 
Theorem 14.11 



Finally, it is sufficient to note that implies j ^{Ci^Ai < 2Rcf3^R J 7(0d?. 
which proves Corollarv 14.21 

5 Conclusions 

It has been shown that a trapped mode is possible in the model presented here. 
To increase confidence that such modes exist on real coasts further work is 
clearly required to demonstrate that this mode is not an artifact of the mod- 
elling assumptions. However these assumptions are the usual ones for the simple 
theory of CSWs and extensions to include stratification and more realistic bound- 
ary conditions have not in general contradicted them [LBMJ. The result here 
suggests that it would be of interest to compare low frequency velocity records 
in the neighborhood of capes with those on nearby straight coasts to determine 
whether there is indeed enhanced energy at the cape. Both the above endeavors 
are being pursued. 
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